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Abstract 

Light-sensitive modification (ruthenium catalysed) of the Belousov-Zhabotinsky 
medium exhibits various regimes of excitability depending on the levels of illumina- 
tion. For certain values of illumination the medium switches to a sub-excitable mode. 
An asymmetric perturbation of the medium leads to formation of a travelling local- 
ized excitation, a wave-fragment which moves along a predetermined trajectory, ide- 
ally preserving its shape and velocity. To implement collision-based computing with 
such wave- fragments we represent values of Boolean variables in presence/absence 
of a wave-fragment at specific sites of medium. When two wave-fragments collide 
they either annihilate, or form new wave-fragments. The trajectories of the wave- 
fragments after the collision represent a result of the computation, e.g. a simple 
logical gate. Wave-fragments in the sub-excitable medium are famously difficult to 
control. Therefore, we adopted a hybrid procedure in order to construct collision- 
based logical gates: we used channels, defined by lower levels illumination to subtly 
tune the shape of a propagating wave-fragment and allow the wave-fragments to 
collide at the junctions between channels. Using this methodology we were able to 
implement both in theoretical models (using the Oregonator) and in experiment 
two interaction-based logical gates and assemble the gates into a basic one-bit bi- 
nary adder. We present the first ever experimental approach towards constructing 
arithmetical circuits in spatially-extended excitable chemical systems. 

Keywords: Belousov-Zhabotinsky reaction, logical gate, adder, unconventional com- 
puter, chemical computer 



1 Introduction 

In 1994 the Showalter Laboratory presented the first ever experimental im- 
plementation of logical gates in the spatially-extended Belousov-Zhabotinsky 
(BZ) system based on the space-time interaction of travelling excitation 
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waves. The logical gates were based on the geometrical configuration of chan- 
nels in which the excitation waves propagated. The ratio between the channel 
diameter and the critical nucleation radii of the excitable media allowed var- 
ious logical schemes to be realised. The original findings by Showalter and 
colleagues led to several innovative designs of computational devices, based 
on geometrically-constrained excitable substrates. 

Designs incorporating assemblies of channels for excitation wave propaga- 
tion were used to implement logical gates for Boolean and multiple-valued 
logic pHl22l)19|r28] . many-input logical gates [T7fT8] . counters [15], coincidence 
detectors [H], and detectors of direction and distance [T61I29] . All these chemi- 
cal computing devices were realised in geometrically-constrained media where 
excitation waves propagate along defined catalyst loaded channels or tubes 
filled with the BZ reagents. The waves perform computation by interacting at 
the junctions between the channels. 

In all papers cited above an excitable BZ medium was employed. In an ex- 
citable medium any initial perturbation will lead to the formation of omni- 
directional wave-fronts of excitation, usually target or spiral waves. In BZ- 
based logical circuits movement of signal, represented by an excitation wave- 
front is controlled by physically restricting medium to dedicated channels of 
propagation. In [2]Pf|27llll|l3] we developed an alternative approach by employ- 
ing self-localised excitations (remotely resembling dissipative solitons) in the 
sub-excitable BZ medium to implement collision-based computing circuits [TJ . 

The paradigm of collision-based computing originates from the computational 
universality of the Game of Life [10] , conservative logic and the billiard-ball 
model [13] with its cellular-automaton implementation [21]. A collision-based 
computer employs mobile self-localized excitations to represent quanta of in- 
formation in active non-linear media. Information values, e.g. truth values of 
logical variables, are given by either the absence or presence of the localizations 
or by other parameters such as direction or velocity. The localizations travel in 
space and collide with each other. The results of the collisions are interpreted 
as computation. There are no predetermined stationary wires, a trajectory of 
the travelling localization is a momentary wire. Almost any part of the reactor 
space can be used as a wire. Localizations can collide anywhere within this 
space. The localizations undergo transformations, form bound states, annihi- 
late or fuse. Information values of localizations are transformed as a result of 
these collisions pQ. 

To implement a collision-based scheme in a spatially-extended chemical medium 
we must employ travelling localised excitations. Such localisations, or wave- 
fragments, emerge in a light-sensitive BZ medium when it is in a sub-excitable 
state [23] • The ruthenium-catalyzed BZ medium shows a high degree of light- 
sensitivity. At some levels of illumination the medium behaves as a classi- 
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cal excitable medium where a perturbation leads to the formation of omni- 
directional propagating waves of excitation. When the level of illumination 
exceeds a critical threshold no excitation persists. There is narrow interval of 
illumination parameters where the BZ medium is in a sub-excitable (weakly 
excitable) state between the non-excitable and excitable states. A perturba- 
tion of the sub-excitable medium leads to formation of localised travelling 
excitations, or wave-fragments. The wave-fragments travel along their prede- 
termined trajectories and preserve their shapes and velocity vectors for some 
period of time. A behaviour of each fragment is determined by the exact level 
of illumination and the size of the fragment. The smaller fragments usually 
collapse, whereas large fragments usually expand. If the illumination level is at 
a critical level then appropriate sized fragments entering a sub-excitable media 
will preserve their size for appreciable distances/time intervals. See overviews 
of recent results in collision-based computing in BZ in piofo] . 

In the present paper we combine the geometrically-constrained approach with 
a collision-based paradigm to deal with the instability of excitation wave- 
fragments. We allow the wave-fragments to travel in channels, defined by low 
levels of the medium's illumination, but collide in larger areas at the junctions 
between the channels. Thus the size of a wave-fragment is "automatically" 
maintained whilst travelling along the channel and enters the junction as a 
localised excitation wave with the ability to expand/contract dependent on 
the level of illumination. The localised excitation waves collide and interact in 
the junctions according to the principles of collision-based computing. 

The paper is structured as follows. In Sect. [2] we present details of the two- 
variable Oregonator model used in simulation of BZ-based logical gates. Func- 
tionality, light-control and polymorphism of simulated ballistic gates are dis- 
cussed in Sect.[3j The Computational model of a one-bit half-adder is presented 
in Sect. |4j In Sect. [5] we outline the experimental techniques used for imple- 
menting BZ-based logical gates. The results of the experimental implementa- 
tion of the interaction gates and one-bit half-adder are discussed in Sect. |6j 
The drawbacks of the current approach and further designs are discussed in 
Sect. 



2 Modelling technique 



We employ the two- variable Oregonator equation [12] adapted for light sensi- 
tivity as an analogue of the Belousov-Zhabotinsky (BZ) reaction with applied 
illumination [9]: 
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The variables u and v represent the local concentrations of activator, or ex- 
citatory component, and inhibitor, or refractory component respectively. Pa- 
rameter e sets up a ratio of time scale for the variables u and v, q is a scaling 
parameter dependent on the rates of activation/propagation and inhibition, / 
is a stoichiometric factor. Constant <fi is the rate of inhibitor production. In 
the light-sensitive BZ represents the rate of inhibitor production which is 
proportional to the intensity of illumination. We integrate the system using 
the Euler method with five-node Laplace operator, time step At = 0.005 and 
grid point spacing Ax = 0.25, e = 0.022, / = 1.4, q = 0.002. The equations 
effectively map the space-time dynamics of excitation in the BZ medium and 
have proved to be an invaluable tool for studying the dynamics of collisions 
between travelling localized excitations in our previous work [2ll4l27|lll] . 



3 Ballistic gates 



Ballistic gates were originally designed to implement collision-based comput- 
ing primitives with propagating slime mould Physarum polycephalum) [?]. 
The idea behind the gate is that if a travelling localisation (slime mould in [7j 
or excitation wave-fragment in present paper) propagates by itself and does 
not interact with other localisations then the localisation always continues 
along its original trajectory without changing its velocity vector. If two local- 
isations collide their velocity vectors are altered and they exit the gates via 
different channels. If the first localisation hits the refractory tail of a second 
localisation then the first localisation annihilates. 

The geometric structure of gates Pi and P2 are shown in Fig. [TJ Input variables 
are x and y and outputs are p and q. Presence of a wave-fragment in a given 
channel indicates Truth, or value '1', and absence — False, or value '0'. Each 
gate implements a transformation from (x,y) — > (p,q). To represent value '1' 
in input channel x or channel y we generate an excitation near the entrance of 
the channel. Two wave-fragments are formed. One of the wave fragment travels 
away from the gate, but the other travels through the channels towards the 
outputs. The fragment travelling within the channels is only considered to 
have registered an output when it reaches the end of the respective output 
channel, we consider the output value in that channel is '1', otherwise if it 
does not reach the end of the channel we consider the output to be '0'. 
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Fig. 1. Geometric structure of gates (a) Pi and (b) P2: x and y are inputs, p and q 
are outputs. From [7j. 

To test the full functionality of the gates, we run simulations with every pos- 
sible input combination (x,y) = (0, 1), (1,0), and (1, 1) as well as a range of 
selected values of G A U B, where: 

• A = {0.07 + n x 0.001 : n is an integer; and < n < 20} 

• B = {0.077 + n x 0.0001 : n is an integer; and < n < 20} 

With the value of <j) increasing from the lowest value used, = 0.07 to the 
largest value, = 0.09, the excitability of the simulated BZ medium changes 
from excitable to sub-excitable, and finally non-excitable. Therefore with the 
same input values, different outputs may be generated. The following is a 
summary of the simulation results obtained when using gates Pi and P%. 



3.1 Gate P x 

For gate Pi, when (x,y) = (1,0), there are four scenarios observed depending 
on the value of 0, as shown in Fig. [2] In the first scenario (Fig. [2^,), when 
< 0.074, the wave started in channel x travels to the intersection and enters 
all three other channels including the other input channel y. This is because at 
this value of the reaction exhibits a high level of excitability. In the second 
scenario (Fig. [2|d), when 0.075 < < 0.0771, the wave reaches the intersection 
of the junctions but only travels to the two output channels p and q. For both 
of these scenarios the output is (p, q) = (1,1). However, you could imagine 
a scenario where this gate with input (x) only would have outputs (p, q, r) 
where output channel (r) is equivalent to input channel (y). 
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(a) (b) (c) (d) 

Fig. 3. Four scenarios observed using gate P\ when input values (x, y) = (0, 1). The 
examples shown here are achieved using specific values of <p: (a) 0.07, (b) 0.076, (c) 
0.0785 and (d) 0.079. 

In the third scenario (Fig. [2^), when 0.0773 < < 0.0777, with the excitability 
of the simulated BZ medium further decreased, only channel q would see an 
output wave-fragment thus the outputs (p, q) = (0, 1). And finally (Fig. |2]l), 
when > 0.0779, the wave-fragment cannot travel far from its origin, due 
to a low level of excitability and also the constraining nature of the channels 
structure, and the outputs are {p, q) = (0,0). 

Similarly, when (x, y) = (0, 1), there are also four different scenarios observed 
in the simulation experiments depending on the different levels of excitability, 
as show in Fig. [3j When < 0.074, the wave-fragment originating from channel 
y travels into the intersection between the channels and into the three other 
channels (Fig. [3^i); when 0.075 < < 0.0772, the wave fragment travels into 
the intersection but only travels to the end of the two output channels p and 
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(a) (b) (c) (d) 

Fig. 4. Three scenarios observed when using gate P\ when the input values 
(x,y) = (1,1) (a-c), and (d) shows an example of an indecisive result. The ex- 
amples shown here are achieved from specific values of 4>: (a) 0.076, (b) 0.0785, (c) 
0.079 and (d) 0.0773. 

q (Fig. (3)3) . Again in both cases the outputs (p,q) = (1,1). Whereas in the 
other two scenarios, when 0.0776 < 4> < 0.0787, the outputs (p, q) = (0, 1) 
(Fig. and when <\> > 0.0788, (p,q) = (0,0) (Fig.gl). 

Finally, when (x,y) = (1,1), only three different scenarios are observed, as 
shown in Fig. gj When <p < 0.0771, (p, q) = (1, 1) (Fig. gi), when 0.0774 < 
<p < 0.0787 (p,q) = (0,1) (Fig. Q)), and when <j> > 0.0788, (p,q) = (0,0) 
(Fig. 

Not every result of the specified values of (ft are included in the above summary. 
For example, as shown in Fig. [4|d, when (x,y) = (1,1) and <p — 0.0773, the 
wave-fragment reached the end of the output channel p, but only part of the 
wave is present due to a lower level of excitability. In this case, it is arguable 
that the output value in the channel p is '1' or '0', and we consider it as 
indeterminate and note it with Thus the output in this particular case is 

(?,?} = (- 0). 

By combining the above results for all possible input values together, we are 
able to determine the functionality of the gate at different values of 0, as shown 
in Tab. [TJ In channel p, one can implement the x V y logic, and in channel q, 
one can implement both the x V y logic and the y logic. 

Furthermore, on closer inspection to Tab. [TJ one may find that for the input 
values of (x,y) = (0,1) and (x,y) = (1,1), the output values are almost 
identical in every range of <fi. This is what we call a y-dominant situation. In 
addition, there is an x-dominant situation. The reason why it is x-dominant or 
^/-dominant, may be seen by studying the time lapse images shown in Fig. |5j 
Both images are achieved at = 0.0779, the difference is that in Fig. [5]d 
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Table 1 

Functionality of gate Pi at different value ranges of <j) (y-dominant). 




(a) 



Fig. 5. Gate Pi has different functionality depending on whether: (a) y-dominant 
and (b) x-dominant. This is determined by which wave- fragment (originating from 
channel x or y) arrives at the junction first. 
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Table 2 

Functionality of gate Pi at different value ranges of <p (x-dominant). 

the origin of the wave in input channel x is farther away from the edge of 
the channel than in Fig. [5^,. And when the two wave-fragments collide and 
join together, it looks as if the wave-fragment from channel x is absorbed, 
thus y-dominant, in Fig. |5^l. And in Fig. [5]d, it looks as if the wave-fragment 
from channel y is absorbed, thus x-dominant. In this case, whether the gate 
Pi appears to be x-dominant or y-dominant is determined by which wave- 
fragment (originating from channel x or y) arrives at the junction first. 

It is evident that in an x-dominant situation, the functionality of gate Pi 
(summarised in Tab. [2]) is different from that when there is a y-dominant 
situation. It shows that for the x-dominant situation, one can implement a 
-ix A y logic. 

3.2 Gate P 2 

For gate P 2 , there were three different scenarios observed. The first one is when 
<p < 0.075 and the simulated BZ medium is excitable, at the junction wave- 
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Fig. 6. Functionality of Gate P2 when the simulated Bz medium is excitable. 
Wave-fragments would arrive at both output channels p and q, thus (p, q) = (1, 1): 
(a) (x,y) = (0,1), (b) (x,y) = (1,0) and (c) (x,y) = (1,1). In all three examples, 
4> = 0.070. 

fragments travel to all possible routes and arrive in both output channels p 
and q in addition to free input channel, as shown in Fig. [6j Note that in Fig. ^jp, 
the wave-fragment originating from channel x arrives at the junction first and 
travels into channel y, where it collides with the wave-fragment originating 
from channel y and the wave fragment originating from y is annihilated. 

In the second scenario, when 0.076 < <p < 0.0787 and the simulated BZ 
medium is sub-excitable, the wave-fragment travels straight through the junc- 
tion to the opposite output channel when there is only one wave-fragment, as 
shown in Fig. [7^, and Fig. [7|d. When (x,y) = (1, 1), when the wave-fragment 
originating from channel y arrives at the junction, it would collide with the 
tail of the wave-fragment originating from channel x and subsequently be an- 
nihilated, as shown in Fig. ^jp. This is due to the relative lengths of the input 
channels prior to reaching the junction. 

In the third scenario, when <ft > 0.0788 and the simulated BZ medium is 
unexcitable, no wave-fragments are observed to arrive at the output channels 
p and q since they collapse not long after being initiated. Therefore the outputs 



(p,q) = (0,0). 

The functionality of Gate P2 with different levels of excitability is summarised 
in Table [U 



4 Simulated one-bit half-adder 

A one-bit half-adder is a logical circuit which takes two inputs x and y and 
produces two outputs: sum xy + xy and carry xy. In previous work a one-bit 
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Fig. 7. Functionality of gate P2 when the simulated medium is close to a sub-ex- 
citable state. Wave-fragments are observed to travel straight through the junction 
to the opposite output channel where possible. The outputs (p, q) are either (1,0) 
or (0,1): (a) (x,y) = (0,1), (b) (x,y) = (1,0) and (c) (x,y) = (1,1). In all three 
examples, 4> = 0.0785. 
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Fig. 8. Functionality of gate P2 when the simulated BZ medium is in an unex- 
citable state. Wave-fragments collapse prior to interacting at the junction therefore 
no excitation waves reach the output channels p and q, thus (p, q) = (0,0): (a) 
(x,y) = (0,1), (b) (x,y) = (1,0) and (c) (x,y) = (1,1). In all three examples, 
6 = 0.079. 



Table 3 



(x,y) 





< 0.075 


[0.076,0.0787] 


> 0.0788 


(0,0) 




(0,0) 


(0,0) 


(0,0) 


(0,1) 


(p,q) 


(1,1) 


(0,1) 


(0,0) 


(1,0) 


(1,1) 


(1,0) 


(0,0) 


(1,1) 




(1,1) 


(1,0) 


(0,0) 


channel p 




X 





channel q 


xVy 


-ix A y 






11 




x< 



xy + xy 



Fig. 9. Scheme of one-bit half-adder constructed from gates Pi and P%. Inputs are 
indicated by arrows. Outputs xy + xy and xy are sum and carry values computed 
by the adder. Outputs and x + y are byproducts. 

half-adder was constructed with Physarum and required two copies of gate P\ 
(Fig. [IJi) and two copies of gate P 2 (Fig. [T]d). Cascading the gates into the 
adder is shown in Fig. [9j Signals x and y are inputted using the two P 2 gates. 
The outputs of the two P 2 gates are connected to the inputs of the two P± 
gates. 



Simulation of the adder using the Oregonator model modified for light sen- 



sitivity was successful (Fig. 10). To simulate inputs x = and y = 1 we 
initiate wave-fragments at the beginning of the y input channels (marked y 
and arrow in Fig. |9|). The wave-fragments propagate along these channels. The 



waves do not branch at the junctions with other channels because we keep the 



wave- fragments localized by varying the parameter (Fig. 10 ib). 



For input values x = 1 and y = wave-fragments are initiated at sites marked 
x and arrow in Fig. |9j The wave-fragment initiated in the left x-input channel 
propagates towards the x+j/-output of the adder. The wave-fragment initiated 



in the right x-input channel travels towards xy + xy (Fig. 10). 



When both inputs are activated, x = 1 and y = 1, the wave-fragment orig- 
inating from the left y-input channel is blocked by the refractory tail of the 
wave-fragment originating from the left x-input channels. The wave-fragment 
traveling in the right x-input channel is blocked by the refractory tail of the 
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(a) 



(b) 




(c) 

Fig. 10. Time lapsed images of localized wave-fragments traveling in the one-bit 
half-adder scheme in the simulated BZ medium. Dynamics of excitation is shown 
for the following input values and parameters: (a) (x, y) = (0, 1) , 4> = 0.0785, 
(b) (x, y) = (1, 0} , <j, = 0.0774, (c) (x, y) = (1, 1) , = 0.0770 



13 



wave-fragment traveling in the right y-input channel. The wave-fragments rep- 
resenting x = 1 and y = 1 enter gate Pi at the top-right-hand-side of the 
scheme and emerge at its outputs xy and x + y (Fig. 10). Thus the functional- 
ity of the designed circuit Fig. 10 is demonstrated in computer simulation. It 
is interesting to note that there is a specific value of required in order to ob- 
serve the correct function of the adder for each possible input sequence. This 
highlights the high degree of controllability that can be implemented when 
using computer simulations. It also highlights in line with the results obtained 
for the single gates how sensitively dependent the functionality of the gate is 
on the level of excitability conferred by <fi. If a single value of was used then 
certain functionality could not be achieved. 



5 Experimental techniques 



Sodium bromate, sodium bromide, malonic acid, sulphuric acid, tris(bipyridyl) 
ruthenium(II) chloride, 27% sodium silicate solution stabilized in 4.9 M sodium 
hydroxide were purchased from (Sigma-Aldrich, U.K., BH12 4QH) and used 
as received unless stated otherwise. 

To create the gels a stock solution of the sodium silicate was prepared by 
mixing 222ml of the purchased sodium silicate solution with 57 ml of 2 M 
sulphuric acid and 180 ml of deionised water. Ru(bpy) 3 S04 was recrystalised 
from the choloride salt with sulphuric acid. Solutions for making the gels were 
prepared by mixing 2.5 ml of the acidified silicate solution and 0.6ml of 0.025M 
Ru(bpy) 3 S0 4 with 0.7 ml of 1M sulphuric acid solution. Using capillary action, 
portions of this solution were quickly transferred into a custom-designed 25cm 
long 0.3 mm deep Perspex mould covered with microscope slides. The solution 
was left for 3 hours to permit complete gellation. After gellation the adherence 
to the Perspex mould is negligible leaving a thin layer on the glass slide. After 
3 hours the slides were carefully removed from the mould and the gels on the 
slides were washed in deionised water at least five times. The gels were 26 mm 
by 26mm, with a wet thickness of approximately 300/im. The gels were stored 
under water and rinsed with deionised water just before use. 

The catalyst free reaction mixture was freshly prepared in a 30 ml continuously- 
fed stirred tank reactor (CSTR), which involved the in situ synthesis of stoi- 
chiometric bromomalonic acid from malonic acid and bromine generated from 
the partial reduction of sodium bromate. This CSTR in turn continuously 
fed a thermostatted open reactor with fresh catalyst-free BZ solution in order 
to maintain a nonequlibrium state. The final composition of the catalyst-free 
reaction solution in the reactor was: 0.42 M sodium bromate, 0.19 M malonic 
acid, 0.52 M sulphuric acid and 0.11 M bromide. The residence time was 30 
minutes. 
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Fig. 11. A block diagram of the experimental setup where A: computer, B: projector, 
C: mirror, D: microscope slide with the catalyst-laden gel, E: thermostatted Petri 
dish, F: CSTR, Gl and G2: pumps, H: stock solutions, I: camera, J: effluent flow, 
K: thermostatted water bath. 



The open reactor was surrounded by a water jacket thermostatted at 20°C. 
Peristaltic pumps (Watson Marlow Ltd. UK. TR11 4RU) were used to pump 
the reaction solution into the reactor and remove the effluent. A diagrammatic 



representation of the experimental setup is shown in Fig. 11 



A Sanyo PROxtrax multiverse projector (Sanyo UK WD24 4PT) was used 
to illuminate the computer-controlled image. Images were captured using a 
Lumenera infinity 2 USB 2.0 scientific digital camera fitted with a macro 
video zoom lens (18-108 f/2.5) and a blue filter (Edmunds optics, UK. Y026 
6BL). The computer generated images of the gates were projected onto the 
catalyst loaded gel via an optically flat mirror (Edmunds optics). The spa- 
tially distributed excitable field on the surface of the gel was achieved by the 
projection of dark coloured channels on a light coloured background. The light 
intensity of the channels was controlled via in house written software. The size 
of the projected channel width was approximately 1.7 mm. Every 10 seconds, 
the pattern was replaced with a uniform light level of 5.7mW cm -2 for 10 ms 
during which time an image of the BZ fragments on the gel was captured. 
The purpose of removing the projected image during this period was to allow 
activity on the gel to be more visible to the camera. 

Captured images were processed to identify chemical activity. This was done 
by differencing successive images to create a black and white image. The levels 
of excitation in the analysed cells are shown in white and the background in 
black. The images were cropped and the channel boundary was superimposed 
on the images to aid analysis of the results. 
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Fig. 12. Structure of the one-bit half-adder showing the position of the interaction 
gate Pi within the overall scheme. 

6 Experimental results 



The scheme for the half bit adder had been successfully designed and imple- 
mented using the Oregonator model of the light sensitive BZ reaction. The 
half adder scheme consists of two pairs of interaction gates (Fig. [I£i,b) which 
must be solved in order to successfully implement the half-adder. Therefore, 
because of the complex nature of the interactions and their high dependency 
on synchronised inputs (more challenging in experimental implementation) a 
strategy of splitting the one-bit half-adder into its constituent parts and solv- 
ing each part separately was devised. This involved projecting only parts of 
the adder scheme onto the thin film light sensitive BZ reaction and injecting 
fragments into the relevant inputs. By varying the light levels in the channels 
the full range of behaviour of the interaction gates could be explored. This 
enabled the optimal light level for solving each gate to be identified. Then 
using the optimal conditions an attempt to construct the one-bit half-adder 
scheme in experiment was undertaken. 



6. 1 Results of implementing interaction Gate P\ 



The geometric structure of gate Pi where x and y are inputs channels and p 
and q are output channels is shown in Fig. [T^i. Also shown in Figure 12 is the 
adder scheme identifying the part of the overall adder which constitutes Pi. 
The structures were drawn in mirror image compared to the gate structures 
shown for the simulation results because in experiment the projected image 
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(a) 

Fig. 13. Progression of fragments seen every 30 seconds when x = 1 at various light 
intensities where (a) excitable at 2.1 mW cm -2 (b) weakly excitable at 2.6 mW 
cm~ 2 (c) sub excitable at 3.1mW cm~ 2 . 

was reflected on to the gel using a mirror. 

In the chemical experiments using the light sensitive BZ reaction and project- 
ing the image of gate Pi, when (x, y) = (1, 0) there are three different types of 
behaviour observed for the interaction of the input fragment at the junctions 
that are sensitively dependent on the excitability level of the channels. The 
excitability level of the channels is directly correlated to the projected light 
intensity. For correct implementation of the half adder scheme where x = 1 



there should only be a fragment travelling into output channel q. In Fig. [13p ,, 
the light intensity was 2.1mW cm -2 meaning that the channels were excitable 
and therefore, the input fragment from channel x is able to enter both output 
channels (p and q) but in addition can also travel into the other input channel 



{%/). In Fig. 13 d the light level was increased to 2.6mW cm -2 making the chan- 
nels weakly excitable. The input fragment travelled into the output q channel 
and input y channel only. This is in contrast to the results observed during 
the simulation experiments where the fragment at an intermediate light level 
{(f) value 0.076 in the model) travelled only to the output channels p and q. 
This difference in observed behaviour is due to a decrease in stability of the 
fragment in the experiment with increasing distance from the initiation site. 
Therefore, at the first junction the fragment expands and is able to transfer 
into the free input channel. However, when the fragment reaches the second 
junction it has lost stability and become localised. The fragment can be seen 
to follow the contour of the straight output channel. 



In Fig. 13:, the light intensity was increased to 3.1mW cm -2 meaning the 
channels were close to the sub-excitable limit. At this light level the fragment 
from the input x channel only travelled into the output channel q at the junc- 
tion between p and q. At this light level the fragments are very weak, diffuse 
and relatively difficult to observe with the camera. Therefore, the differenced 
image does not show the fragment travelling through the entire output chan- 
nel. 

When (x, y) = (0, 1) there were also three distinct types of behaviour observed 
at different light levels. The desired behaviour in order to implement the adder 
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Fig. 14. Progression of fragments seen every 30 seconds when y = 1 at various light 
intensities where (a) excitable at 0.4mW cm -2 (b) weakly excitable at 2.6mW cm~ 2 
(c) sub excitable at 3.1mW cm~ 2 . 




Inpul Input 



Fig. 15. Experimental results when X = 1 and Y = 1 at 2 different light intensities 
(a) 2.1mW cm~ 2 and(b) 3.1mW cm~ 2 . 

successfully was that an input fragment in the y channel only would result in 



a single fragment in the q output channel. In Fig. 14 1 the light intensity was 
0.4mW cm -2 in the excitable domain and again this resulted in the input 
fragment travelling in the y channel splitting into both output channels p and 



q and also into the input channel x. In Fig. [14p the light level was increased 
to 2.6mW cm -2 in a weakly excitable domain and the fragment travelling 
in the y input channel split at the junction and travelled into the q output 
channel and the x input channel. Again this is in contrast to the simulation 
results where the fragments travelled to the output p and q channels only at 



the intermediate light level. Finally in Fig. [14p the light intensity was 3.1mW 
cm -2 the channels were sub excitable and the fragment travelled only into the 
output q channel. 

When (x,y) = (1, 1) the desired behaviour (in order to solve the half adder) 
is for the fragments to collide and merge at the intersection of the output 
channels and then split to give a fragment in both the p and q output channels. 
At excitable light levels up to a value of 2.1mW cm -2 the fragments travelled 
from the input x and y channels collided, merged and split at the junction to 



give fragments at the output p and q channels (Fig. 15 %). However, if the light 
level was increased to 3.1mW cm~ 2 close to the sub-excitable domain then 
despite collision of the input fragments a fragment was only observed to travel 
into the q output channel (Fig. [15p). Therefore, despite the optimal light level 
for solving the (x,y) = (1,0) and (x,y) = (0,1) cases being 3.1mWcm~ 2 
this is different to the light level found to be optimal for solving the (x, y) = 
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Fig. 16. Adder structure showing the position of interaction gate Pi- 

(1, 1) case. This may have implications for solving the adder using the current 
scheme if the light level is held at a constant value. 

Despite the fact that this interaction gate is part of a larger adder scheme 
the results obtained in experiment for the (x, y) = (1,0) or (0,1) highlight 
the adaptability of using light to impart different functionality to logic gates. 
For both the x = 1 and y = 1 case the same general behaviour is observed. 
Low light intensity (excitable), x = 1 or y = 1 output 1, 1 and 1. Medium 
light intensity (weakly excitable), x — 1 or y — 1, output 1,1 and 0. High 
light intensity (sub-excitable), x — 1 or y — 1, output 1,0 and 0. Although the 
fourth case is not investigated it is obvious that if the light level is increased 
further then the x — 1 or y — 1, output 0,0 and could be obtained in 
experiment. Also in our experimental results we count any excitation at the 
end of the output channel as 1. However, as demonstrated for the simulation 
experiments there can also be an indeterminate result. This could be achieved 
in experiment by setting a threshold range for the level of excitation detected 
at the output channel. In this way the functionality of the gate could be further 
increased at differing light levels. In experiment it is not as simply to map all 
possible light levels because the system of changing the light level does not 
have as many degrees of freedom as changing parameter in the simulation 
experiments. 

6.2 Experimental results for interaction Gate Pi 

The geometric structure of gate P2 where x and y are input channels and p 
and q are output channels is shown in Fig. [TJd. The adder structure showing 
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Fig. 17. Progression of fragment every 30 seconds when x = 1 at the cross junction 
(a) excitable at 2.2mW cm -2 (b) sub-excitable 3.1mW cm -2 . 



q 



X 





Fig. 18. Geometric structure of gate (a) P 2 , where x and y are inputs, p and q are 
outputs, (b) Adder structure showing the position of gate P 2 . 



the position of gate P2 is shown in Fig. [16} During these experiments only a 
single section of the adder scheme is solved where the y input channel is longer 
than x. In the other half of the adder scheme this is reversed and the x input 
channel is longer than y. For gate P 2 to be solved correctly an optimal light 
level has to be found so that a fragment in either the x or y channel will travel 
straight into either the q or p output channels respectively. When there are 
inputs at both x and y the dominant fragment (from the short input channel, 
in this case x) should block the fragment from the longer input channel (in 
this case y). The resulting collision at the junction should result in a fragment 
travelling into the p output channel (in this case). 



For gate P 2 it was found in experiment that at a lower light intensity of 
2.2mW cm -2 for the (x, y) = (1, 0) case the fragment from the input x travels 



to the input y and both output channels p and q see Fig. [T7p ,. This is because 
at this light level the channels are in an excitable domain. However, if the 
light intensity is raised to 3.1mW cm -2 close to the sub-excitable limit then 
the fragment from the input x channel only travels to the output p channel 
(Fig.[fi}>). 
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Fig. 19. Progression of fragments every 30 seconds when y = 1 at the cross junction 
at (a) 3.7mW cm" 2 and (b) 4.4mW cm -2 . 




Fig. 20. Progression of fragment every 10 seconds when x = 1 and y = 1 at light 
intensity 1.5mW cm -2 (excitable) (a) x fragment is dominant (b) this results in a 
collision and the annihilation of the y fragment (c) the resulting daughter fragment 
spreads across the north and east channels (d) finally it breaks into 2 fragments 
which travel along the north and east channels (grey arrows show direction of frag- 
ment input). 



For the case (x,y) = (0, 1) the interaction gate P' 2 was used (Fig. [18ft), note 
also the position in the overall adder scheme in Fig. [T8| 3. At a lower light 
intensity of 3.7mW cm -2 close to the sub-excitable limit the fragments from 
the input y channel travels into the input x channel and both output channels 
p and q (Fig. 19 1). If the light level is raised further to 4.4 mW cm -2 then 
the fragment only travels into the q output channel (Fig. [lOpj). The light 
levels are different from the previous experimental results for (x, y) = (1,0) 



(Fig. |17p,b) due to subtle differences in gel composition and morphology which 
can make the sub-excitable region variable in experimental implementations. 
However, the sub-excitable region can always be identified by the formation of 
unbounded fragments in a projected uniform light level prior to the projection 
of the interaction gates/channels. 



Figure 20 i-d shows the results for (x, y) = (1, 1) at a light intensity of 1.5mW 
" 2 in the excitable domain. In this implementation the fragment x is dom- 



cm 



inant (despite the y input channel being shorter) and it undergoes collision 
with fragment y, resulting in annihilation of the y fragment. However, due to 
the excitability of the channels the fragment expands in a north east direction 
at the junction and splits into two fragments that travel along the output 



p and q channels. Figure 21i-d shows the converse case where the y frag- 
ment is dominant and annihilation of the x fragment occurs and the resulting 
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Fig. 21. Progression of fragment every 10 seconds when x = 1 and y = 1 at light 
intensity 3.0mW cm~ 2 (weakly excitable) (a) y fragment is dominant (b) so anni- 
hilates the x fragment as they collide (c) then spreads across the north and east 
channels (d) then breaks into 2 fragments which travel along the north and east 
channels (grey arrows show the direction of the fragment input). 





(a) 

Fig. 22. When x = 1 and y = 1 at light intensity 3.1mW cm -2 (a) both fragments 
collide to give a V-shaped fragment travelling in a north east direction, progression 
of fragment every 30 seconds (b) x fragment is annihilated by the y fragment and the 
resulting daughter fragment travels north into the q output channel only, progression 
of fragment every 20 seconds. 

daughter fragment also travels towards north east direction then breaks into 
2 fragments which travel along the output p and q channels at light intensity 
3.0mW cm -2 (channels are weakly excitable). 



Figure 22 1 shows the case when both input fragments collide at the same 
time at the cross junction, collide to form a daughter fragment which splits 
into two fragments at the cross junction and travels along the output p and 
q channels even when the excitability is relatively low (light intensity 3.1mW 
cm -2 ) (Fig. 22 1). However, at the same light level if the y fragment is slightly 



dominant then the x fragment is annihilated and the resulting daughter frag- 
ment travels into the output q channel only, see Fig. 
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This is the correct 

implementation for the half adder scheme to be solved. However, this em- 
phasises the importance of synchronising the timing of the two inputs in the 
experiments. This is difficult relative to the simulation experiments although 
differences in the position of the initiations relative to the input channels gave 
similar results in the simulation experiments. In the experiment we utilise a 
constant fragment generator in a dark region which is then linked to two same 
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Fig. 23. Progression of fragments seen every 30 seconds at sub excitable level when 
(a) x = 1 (with one x input only) at light intensity 3.1mW cm -2 and (b) y = 1 
(with one y input only) at light intensity 3.1mW cm -2 , (c) x = 1 and y = 1 (with 
fragments at both x inputs and both y inputs simultaneously) at light intensity 
3.1mW cm~ 2 . 

length excitable channels that transfer fragments to the inputs. These chan- 
nels are removable and in theory should deliver one fragment to each input 
at the same time interval. However, even a small difference in timing has a 
marked effect on the operation of the cross-junction as the lengths of the x 
and y input channels are supposed to confer the functionality of the gate. This 
error in timing is compounded by the fact that the input channels are only 
a few mm in length. The results do show that a wide range of behaviour is 
possible for the same junction by altering the light level and input timing. 
At a critical light level and with the correct timing (i.e. both fragments en- 
ter their respective input channels at the same time) then the cross junction 
is successfully implemented according to the half adder scheme (where they 
must provide the appropriate inputs for gates PI). 



6.3 Attempt to construct one-bit half-adder in experiment 



Figure [23] shows the progression of fragments seen every 30 sec at a sub- 
excitable level, light intensity 3.1mW cm -2 when (a) (x,y) = (1,0) with one 
x input only and (b) (x,y) = (0, 1) with one y input only, (c) (x,y) = (1, 1) 
with fragments at both x inputs and both y inputs simultaneously. As previ- 
ously observed at this light level when solving the interaction ates PI and P2 
an input in either channel x or y results in an output in either channel p or q 
respectively thus giving the correct output according to this theoretical con- 
struct of a one-bit half-adder. Although only one x and y input was initiated 
for ease of experimental implementation each part of the junction has been 
solved previously so solving two simultaneous x or y inputs in experiment at 
this light level should be achievable. For the case where x = 1 and y = 1 both 
sets of inputs were initiated simultaneously in order to implement the half bit 
adder. Both respective cross junctions were solved successfully in this exper- 
imental implementation. However, despite a collision prior to the two output 
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channels a fragment could only be observed in one output channel q. There- 
fore, the entire half bit adder scheme is not solved at a single light level in 
experiment. If we revisit the fractional gate implementation in experiment we 
see that this is actually in agreement with the results from Fig. 15 where the 
collision sequence where x = 1 and y — 1 is only solved at higher excitability 
levels. Therefore, the half bit adder could be solved in experiment via projec- 
tion of differential light levels (or by altering the light level online) in order 
to achieve all parts of the adder scheme in combination. If we compare this 
to the simulation results each part of the adder x = 1, y = 1 and x = 1 and 
y = 1 are solved at subtly different values of all close to the sub-excitable 
limit. In experiment we don't have such fine control over the light dependent 
factor meaning that even if a single light level to solve the adder scheme when 
x and y = 1 does exist (as predicted by simulation) the degrees of freedom in 
altering the light level make finding it improbable especially when considering 
the heterogeneity of the chemical system (which has limited controllability). 



7 Summary and discussion 



We demonstrated in numerical simulations that the light-sensitive BZ medium 
near a sub-excitable level realizes a wide variety of logical operations when 
weakly-constrained by illumination-induced geometrical architectures. We de- 
signed two types of Boolean logic gates, both gates have two inputs and 
two outputs. The gates implement transformations (x,y) — > (xy,x + y) and 
(x,y) — > (x,xy). In addition by mapping all possible light levels and input 
states of these interaction gates a number of intermediate output states were 
identified showing that complex "light mediated" logic gates can be con- 
structed from relatively simple interaction gates. We shown how the inter- 
action gates can be assembled into a one-bit half-adder. The full functionality 
of the adder is illustrated using a two-variable Oregonator model modified 
to account for light sensitivity. The different input states of the adder were 
implemented using different values of0. This fine control of (ft prevents wave- 
fragments from expanding in the propagation channels and junction-chambers 
in order to confer the correct functionality of the integral parts of the adder. 
This is the first instance of arithmetic circuits being implemented in models 
of excitable chemical systems. In laboratory experiments we were able to im- 
plement the interaction gates that make up the one-bit half-adder scheme. We 
were able to show that in agreement with simulation results the functionality 
of the simple gate is sensitively dependent on the light level. Therefore, for the 
same input states we were able to show various output states dependent on 
the light level. Even though some of the functions of the interaction gates dif- 
fered between experiment and simulation they were in close agreement when 
considering the functionality of the interaction gate needed to solve the adder 
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effectively. Therefore, in both the simulation and experimental implementa- 
tions a high level of illumination close to the sub-excitable limit produced a one 
input/one output state for both gates PI and P2. In addition they produced 
a two input /one output state for gate p2 provided the input timing was cor- 
rect. For gate PI the correct implementation was two input/two outputs but 
this was only achieved at relatively low light levels in experiment compared to 
simulation experiments. In laboratory experiments to implement the one-bit 
half-adder there was not a single light level that gave the correct implementa- 
tion. In some respects this agrees with simulation results as for different input 
sequences different levels of are required to stabilise the wave fragments. 
Therefore, there is ot a single value of which gives the correct outputs for all 
possible input sequences. However, in the simulation all the light levels used 
are close to the sub-excitable level whereas in experiment to implement the 
final collision sequence in gate PI requires higher levels of excitability (low 
light levels). Therefore, for the interaction gates PI and P2 to work in com- 
bination the light level would have to be projected differentially or changed 
online in order to implement the adder correctly. Therefore, work will continue 
to identify experimental parameters or different gate schemes whereby arith- 
metic circuits can be implemented effectively in experiment. These laboratory 
experiments have shown that the control of propagating wave-fragments is ex- 
tremely difficult due to the fragments very high sensitivity to environmental 
conditions. In this case the difficulties were multiplied as solving the adder 
scheme required precision timing of the inputs and the sub-excitable nature 
of the chemical media meant that the fragment stability was so low that the 
whole scheme had to be solved at a mm scale. Therefore our future task will be 
to increase the stability of the wave-fragments and to realise correct function- 
ing arithmetic circuits in experiment and many-bit full adders in simulation 
and experiment. 
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